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Abstract 

The CP^ spin model is simulated at large correlation lengths 
in two dimensions. An overrelaxation algorithm is employed 
which yields reduced critical slowing down with dynamical expo- 
nents z ^ 1. We compare our results with recent multigrid data 
on the massgap m and the spin susceptibility and confirm the 
absence of asymptotic scaling. As a new result we find scaling for 
the universal topological susceptibility with values extrapolating 
to Xt/^'^ = 0.156(2) in the continuum limit. 
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1 Introduction 



With the discovery of Monte Carlo simulation algorithms with reduced crit- 
ical slowing down it has become possible to nonperturbatively study a num- 
ber of two dimensional lattice field theories close to the continuum limit. 
As the ratio between physical length scales and the cutoff length reaches 
values of O(IOO), these high precision simulations obtain a completely new 
quality compared with what could be done before. Cluster algorithms ^ 
are particularly efficient, but so far their success has remained restricted to 
0{n) invariant vector models as far as continuous field manifolds are con- 
cerned. Recently multigrid techniques succeeded in extending the class 
of models accessible to such high precision study to include the CP^ and 
SU{3) X SU{3) systems with standard action [If. In parallel with these 
nonlocal methods overrelaxation (OR) techniques have been developed for 
practical applications to spin systems [|, ^ . Although OR is expected to only 
reduce critical slowing down to dynamical exponents 2 ~ 1, thus producing 
a gain in efficiency proportional to the correlation length, it already makes 
a large difference in practice. Being a local algorithm it can be implemented 
very efficiently, with no overhead, and it is actually superior to multigrid for 
many simulations at intermediate correlation lengths of 0(10). In this letter 
we present OR results on the CP^ model. The purpose is on the one hand 
to offer an independent check on results obtained with the new multigrid 
technique in and on the other hand to supplement these data by values 
of the topological susceptibility, a particularly interesting observable in CP" 
modelsQ. In sect. 2 we define the model and present the OR algorithm used. 
The physical results extracted for the CP^ model are discussed in sect. 3. 



2 CP^ spin model and overrelaxation 

The spins of CP^~^ models live on the sites x of a hypercubic L x L lattice and 
consist of rays in complex n dimensional space. They can be labelled by unit 
vectors z{x) entering into the standard action only through 1-dimensional 

^ During our numerical work we received a copy of |^ where |Q is refined by higher 
statistics, and the topological charge is studied, too. The results are now in slight conflict 
with ours, as we discuss in sect. 3 
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projectors 

P{x) = z{x)z\x). (1) 
The partition function then reads 

Z= /nd"'2(a;)dV(a;)(5(|z(a;)|=^ - 1) exp|2/3^tr[P(a;)P(a; + /i)]| , (2) 

where /i — 0,1 and /t are unit vectors in the corresponding lattice directions. 
The fundamental correlation function is given by 

C{x) = (tr[P(x)P(0)]) - 1/n , (3) 

and the magnetic susceptibility x is constructed by summing the argument 
of C over all sites. Wc employ periodic boundary conditions, and the mass 
m or correlation length ^ = are defined through the temporal decay of 
C at zero spatial momentum [we use lattice units putting the lattice spacing 
a = 1]. 

Topological quantities in the CP""'^ model are associated with the "com- 
posite" U{1) gauge field 

_ z^{x)-z{x + fi) 
'~ \z^x)-z{x + ft)\- ^ ^ 

As one forms plaquette fields Up from it and Fp with 

Up = e'^" , -TT < Pp < TT , (5) 

the topological charge 

Q=h^F. (6) 

is integer, since the product over all Up is unity on the torus. While the 
average of Q vanishes for symmetry reasons (which we monitored) , the topo- 
logical susceptibility 

Xt = {Q')/L' (7) 

is expected to be a nontrivial physical quantity of length dimension —2. 

Local updates of a CP"~^ spin z at some site x have to sample the local 
Boltzmann factor 

p{z) oc ex]){z^ M{x)z), (8) 
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where the hermitean n x n matrix M is given in terms of the nearest neigh- 
bors, 

M{x) = 2p J2 Piy)- (9) 

\y-x\ = l 

The decisive component in our reahzation of OR are microcanonical moves 
of individual spins which leave unchanged. As has already been observed 
in 0, this can be achieved, if a normalized eigenvector ip of M{x) is known, 
by reflecting 

-z + 2^^^ ■ z . (10) 

In the CP^ model, which coincides with the 0(3) model, the standard micro- 
canonical reflections P] are reproduced for either choice of the two possible 
eigenvectors of M. For the case n = 4, which we study here, we use the 
eigenvector belonging to the largest eigenvalue, as that one is easy to obtain. 

In the practical realization of hybrid OR we perform microcanoni- 
cal sweeps with checkerboard ordering followed by one standard ergodic 
Metropolis sweep. In the latter, we move the z spins with U{2) rotations 
of all possible pairs of components of z. The eigenvectors of the M{x) are 
simply constructed approximately by multiplying M a few times on a start 
vector preconditioned to the sum of the nearest neighbor z{y) spins. Then 
the reflected spin ( p!0D is first considered as an update proposal, which is 
followed by an accept/reject step taking into account small changes in en- 
ergy due to the approximation. With four multiplications we find for all our 
simulated (3 values between 96% and 99.7% acceptance (growing with /3), 
which makes this step practically microcanonical. Note that in this way the 
remaining small error in the eigenvectors does not lead to any systematic 
bias in the simulation. 

The resulting integrated autocorrelation times as functions of the corre- 
lation length ^ are summarized in Fig.l. Errors are displayed but mostly fall 
within the symbols. We employed the tuning strategy of taking the number 
of microcanonical sweeps N roughly proportional to ^. This was found to 
be optimal both in an extensive numerical study of the 0(3) model and 
in a rigorous analysis of the Gaussian free field model |^ . More precisely, in 
the order of growing ^, we took iV = 1, 3, 6, 9, 12, 25. All Tint refer to single 
sweeps, and we did not distinguish OR and Metropolis here. In fact, close to 
the continuum limit, only the microcanonical OR steps, which are faster, are 
important. The dynamical exponents z shown in Fig.l result from fits to the 
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Figure 1: Integrated autocorrelation times in single sweeps for energy E, spin 
susceptibility x, and topological susceptibility Xt- 

four largest ^ values. While the zero momentum quantities x Xt show 
the expected ^ ~ 1, the value for the energy seems to be even smaller. As E 
can be written with a lattice derivative and couples to all nonzero momenta, 
this behavior cannot be excluded for an integrated autocorrelation time even 
if the slowest modes evolve with z ^ 1. The analogous situation for free 
fields is discussed in P . The topological charge was found to have very short 
autocorrelations, which coincides with observations in P, ^. 

With the results on exponents 2; in 0], it is clear that for asymptotically 
large correlation lengths multigrid becomes superior to OR. Based on CPU 
times kindly communicated by the authors and on our own, we estimate 
the crossover in efficiency to occur at ^ ~ 20-30. In a version of hybrid 
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OR was tested which makes microcanonical moves in the same embedded 
U{1) spin models as are used for multigrid. The crossover occurs at similar 
^, and exponents z close to unity are reached. It is interesting that the 
restriction of microcanonical moves to a small submanifold of spin space 
seems to be of no disadvantage. This opens the possibility of very simple 
exactly microcanonical updates in complex field spaces like SU{3). 

Let us finally mention that the data shown here represent a total of 166 
hours on one XMP processor, of which 2/3 went into the largest lattice. 



3 Results for the CP^ model 

Apart from ^ a number of other studies of the CP^ model have been 
attempted recently, among them |0, |]. They both used cluster algorithms 
that turned out not to be really efficient enough to obtain high precision 
data at large correlation length. The reason for the choice of CP^ was in 
most cases the expectation that pathologies plague lattice definitions of 
the topological charge for smaller n, at least when the standard action and 
lattice are used. 

Results from our runs on lattices are now summarized in the table: 



p 


L 


X 


Xt X 103 


6 


6 


2.3 


20 


11.09(3) 


11.11(4) 


2.579(5) 


2.613(13) 


2.5 


32 


26.17(5) 


5.20(2) 


4.442(7) 


4.488(16) 


2.7 


64 


79.4(2) 


1.731(8) 


8.800(21) 


8.845(46) 


2.8 


96 


145.9(5) 


0.893(5) 


12.640(36) 


12.784(84) 


2.9 


128 


274.7(1.5) 


0.435(4) 


18.49(9) 


18.79(22) 


3.1 


256 


942.6(7.5) 


0.103(1) 


38.37(25) 


39.32(59) 



All la errors quoted have been determined by a jackknife binning proce- 
dure with various bin lengths and typically a few hundred effectively indepen- 
dent bins. This has been compared with the directly summed autocorrelation 
functions for E, x, Xt- The two columns for the correlation length ^1^2 require 
some discussion. We stored all estimates of the two point function at time 
separations t = 0, A,2A, . . . with A = L/32 or A = 1 for L = 20. From 
two successive such values, time dependent effective masses were determined. 



5 



The masses whose time separations selfconsistently embrace are quoted 
as ^j. The table shows that ^i, for our statistical accuracy, is systemati- 
cally shorter than ^2- This is the expected signature of higher mass states in 
the channel probed by @. Inspection of ^3 suggests that, with the present 
level of errors, ^2 is an acceptable estimate for the inverse massgap. The 
ratio ^3/^2 is consistent with unity at the la level, although the mean ex- 
ceeds unity for all but the last line in the table. We face here the recurring 
dilemma of having to compromise between systematic and statistical errors 
when estimating energies. As a consequence, it seems hard to exclude that 
under a worst case scenario, with two nearby levels in the probed channel, we 
still underestimate systematic errors. Comparing with cluster algorithms at 
this point, we see that OR (and multigrid) is inferior not only in speed, but 
also the lack of improved estimators is a severe disadvantage that is costly 
to overcome with just statistics. Finite size effects, on the other hand, are 
expected to be irrelevant with mL ^ 8, and this has been tested to some 
degree in [§. Note that with the above procedure of extracting masses, and 
for the chosen /3 values, all scales (L, ^, A) of the various lattices are roughly 
in the same proportion to one another. Therefore, neglecting scaling viola- 
tions due to the finite lattice spacing, the ratio ^2/^1 should be universal to 
a good approximation. Our data are consistent with this. If we assume it, 
on the other hand, we can determine the ratio to be about 1.01 from the 
smaller lattices, and apply this correction to ^1 on the largest two lattices. 
This gives a smaller error than taking ^2, and we feel that the neglect of 
scaling violations for this 1% correction is no problem. We use these mass 
estimates in the following analysis. 

In asymptotic scaling was found to be violated in a way similar to 
or even worse than for the 0{n) models |jll|. We confirm this and plot in 



Fig. 2 the ratio of the mass and the perturbative lattice scale in the two loop 
approximation 

= (npy/h-^P. (11) 

Asymptotic scaling requires it to tend to a constant as with a speed 

linear in j3~^. Clearly, no sign of this is visible in the (small) range over 
which we are able to vary j3~^. We also consider a modified bare coupling as 



discussed in [11 1. From the perturbative expansion of the energy 



E = {tT[P{x)P{x + A)]) = 1 - - + 0(^3) (12) 

4 Id 
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Figure 2: Ratio of the mass to the two loop approximation of the perturbative 
lattice scale 



we find that 

= 4(1^ (1^' 

is a correctly normalized alternative inverse bare (= short distance) coupling. 
Now A.\ can be reexpanded in 13^ , and the series is truncated to two loops 
in this parameter, whose values are taken from the measured energies. This 
leads to the other data points in Fig.2. They clearly vary less steeply, just like 
in the 0{n) case, but the behavior is not monotonic, and no really convincing 
extrapolation to =0 is feasible either. The difference between the two 
sets of data is best taken as evidence of how far away we are from asymptotic 
behavior in the bare coupling. In P] a different lattice regularization of the 
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Figure 3: Scaling behavior of the topological susceptibility in physical units 

CP^ model has been employed, and asymptotic scaling seems to be violated 
less drastically. It would be very interesting to develop OR or multigrid for 
this action and investigate universality for this case. 

The situation with regard to scaling is rather different for the dimension- 
less product of physical quantities Xt^"^- In Fig- 3 it is plotted against 1/^^, 
the squared lattice spacing in physical units. We see that data from the four 
largest lattices fit with a straight line and extrapolate to 0.156(2). The in- 
trinsic length scale for the scaling violations seems to be about three lattice 
spacings, and for such correlation lengths their expansion in powers of the 
lattice spacing ceases to be useful. The breakaway is rather abrupt. While 
for ^ = 8.8 the leading term alone works very well, each of the smaller lattices 
needs a new power of to be included to get an acceptable fit. We wish 
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to remark that an asymptotic expansion in powers ^ (neglecting the vari 



ation of logarithms) is suggested by renormalized perturbation theory [|T2 
and by an effective action picture [^]. From the point of view of statistical 
mechanics, an arbitrary power of would seem natural. We found that 
the four leftmost points in Fig. 3 can indeed be extrapolated with anything 
between to As our data get very close to the continuum in Fig. 3 
(compare this with Fig. 2!), the extrapolated values are not too different. If 
one wants to include these uncertainties in a systematic error one could still 
quote Xtlfn^ = 0.16(1) in the continuum. 

In xtC,"^ is found to fall again at P = 3.1,3.3 after rising and agreeing 
with our data for smaller f3 values. Since we do not have data for /3 = 3.3, 
the only real discrepancy in the measured data hinges upon a smaller value 
beyond errors for ^ at = 3.1 in 0. As mentioned in the authors 
extract correlation lengths by fitting over a range of separations starting at 
^. We feel that there may also be a problem with higher states and propose 
a simple analysis with effective masses. In any case this discrepancy should 
be resolved in the future, if necessary, with new long runs. 

Let us finally compare with results from the large n expansion [|l^. To 
leading order one expects 

Xt^^ = = 0.05968 ... for n = 4, (14) 

Ann 

from which we differ by a factor of about 2.6. First 1/n corrections have 
been computed in [|14|. The authors find 



1 / 38 \ 

Xti'^ = 1 - — I = 0.03979(1 - 0.095) for n = 4, (15) 

27m V n J 

where ^ is a second moment definition of the correlation length. The authors 
argue that, although = 2^^/3 holds at n = oo, the two lengths differ only 
by a few % for small n with the connection being nonanalytic in 1/n. This 
leads to an even larger discrepancy between our result and large n despite 
an already small leading correction for n = 4. In summary, no quantitative 
agreement with the 1/n expansion at n = 4 is found. 
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